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Abstract 
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1 Introduction 



Monte-Carlo simulations in lattice QCD enable one to calculate the spectrum of low-lying hadrons 
from first principles and thus to obtain important information about the long-range dynamics of quarks 
and gluons. In the past, a major effort has been concentrated on the study of the ground-state particle 
spectrum. However, recently there has been a lot of work on the spectroscopy of the excited nucleon 
states as well [1-9]. For the status of lattice calculations of the baryon spectrum, see e.g. the recent 
reviews [10, 11]. We note that several long-standing puzzles in the field are still awaiting a resolution. 
In particular, it is worth to mention the problem of the level ordering of the negative-parity nucleon 
excitation iV* (1535) and the positive-parity Roper resonance N* (1440) as well as the structure of the 
A (1405). 

The A-resonance is the most important baryon resonance. Its mass is close to the mass of the 
nucleon, and it couples strongly to nucleons, pions and photons. It is clear that a systematic study of 
the properties of the A-resonance in lattice QCD could lay a solid theoretical basis for understanding 
the low-energy QCD dynamics in the one-baryon sector in lattice QCD. Evaluating the mass of the 
A-resonance has already been addressed in the last few years. For illustration, we consider the recent 
calculation of the A-mass in quenched QCD using a tadpole-improved anisotropic action and FLIC 
fermions [8]. In this paper, the local interpolating field with quantum numbers of the A + has been 
chosen in the following manner 

Xtl = ±= eabc j 2 (u aT C 7 ^V + (u aT C 7 ^V } , (1) 

where a, b, c denote color indices and C is the charge conjugation matrix. With this interpolating 
field one further calculates the two-point correlation function at zero spatial momentum and at large 
Euclidean time, extracting the mass spectrum. More precisely, in Ref. [8] the calculation of the ratio 
m^/misr (as well as the similar ratios for some other excited baryon states) has been performed down 
to the quark masses corresponding to (M w /M p ) 2 = 0.4. At such quark masses, the dependence of the 
data on (M^/Mp) 2 is very smooth (almost linear). However, a large gap in the quark (pion) mass has 
still to be bridged if the data are extrapolated down to the physical value of {M^/Mp) 2 ~ 0.03. 

The above example highlights all features of the lattice calculations of the parameters of baryon 
resonances (and the A-resonance, in particular), which we wish to put under scrutiny in the present 
paper: 

a) As already mentioned, the distance from the lowest data point to the physical value of the 
quark mass is still very large and this will lead to large extrapolation errors. In case of a stable 
particle, one would argue that performing calculations at smaller quark masses will reduce this 
uncertainty until, eventually, the simulations are done close to or at the physical value of the 
quark mass. However, in the case of a resonance, the threshold value of the quark mass exists 
after which e.g. the A starts to decay into a pion and a nucleon. It is evident that, below the 
threshold, the method which was described above can not be applied straightforwardly. Does 
this mean that, in the case of a resonance, there is always a (not so small) gap in the lattice 
data, where one should rely only on chiral extrapolation? 

b) Lattice data are always real. Does this mean that one gets the real part of the resonance pole 
mass as a result of a chiral extrapolation below the decay threshold? 

c) Can one determine the decay width of a resonance by combining the method described above 
with the chiral extrapolation? 



2 



Note also that the Monte-Carlo simulations e.g. in Ref. [8] have been carried out in the quenched 
approximation. It is however clear that unstable systems at small quark masses, which will be con- 
sidered below, can be meaningfully discussed only in the context of lattice data based on simulations 
with dynamical fermions. 

To summarize, we want to ask whether there exists a modification of the above method that enables 
one (at least, in principle) to carry out calculations of the properties of an excited state at such quark 
masses when this state becomes unstable. Moreover, one would like to eventually walk all way down 
to the physical quark mass, excluding the extrapolation error altogether (like this appears possible for 
ground-state hadrons). It would also be very instructive to see in detail what happens in the vicinity 
of the threshold, when one crosses it while performing a chiral extrapolation of the spectrum. Finally, 
one should investigate which quantities are most sensitive to the resonance parameters at small quark 
masses, and how accurately one could extract the mass and the width of the resonance by calculating 
these quantities on the lattice. 

The question of identification of hadron resonances on the lattice has already been addressed in 
the past. For example, we would like to mention the papers [12-19]. In general, the method proposed 
originally by Liischer considers the extraction of the two-body scattering phase shifts from the energy 
levels, calculated in a finite Euclidean box. In particular, the signatures of unstable particles have 
been studied. It is demonstrated that in the presence of a narrow resonance the dependence of the 
energy spectrum of the system on the box size L exhibits a very peculiar behavior near the resonance 
energy, where the so-called "avoided level crossing" takes place, see Fig. 1. Note that usually this 
name is used to describe an abrupt rearrangement of the structure of the energy levels of a system, 
which takes place near the resonance energy when the scattering phase passes through it/ 2. From 
Fig. 1 it is clear that the position of a narrow resonance can be readily identified by measuring the 
energy levels of a system at a finite volume and locating horizontal "plateaus." Moreover, the minimal 
distance that separates the energy levels near the avoided crossing, is determined by the decay width 
of a resonance and, consequently, the latter quantity can be also extracted from the same lattice data. 
A very nice qualitative discussion of the avoided level crossing is given e.g. in Refs. [15, 16,20]. 

An approach, related to the Liischer's method, is used in Ref. [20] to study the structure of the 
energy levels in the two-pion system carrying the quantum numbers of the p-meson. To this end, these 
energy levels are calculated in an effective field theory (EFT) with a phenomenological lowest-order 
pirn coupling. Furthermore, in Ref. [21] it has been shown that the presence of a narrow excited state 
above the threshold modifies the simple exponential decay law of the time-sliced two-point function. 
The decay width within this approach is extracted not from the two-point function, but directly from 
the decay amplitudes (see also [22,23]). In addition, we point out the recent investigation [24], where 
it is proposed to reconstruct the spectral density in the two-point function by using the maximum 
entropy method. This approach, in principle, also has the capability to address the problem of unstable 
states in lattice calculations. 

Note also that Liischer's approach has been recently applied to study nucleon-nucleon phase shifts 
at low-energy, as well as two-body shallow bound states [25-29]. To this end, in Ref. [25], Liischer's 
master formula which relates the scattering phase shifts to the energy level displacements caused by 
the NN interaction, is re-derived within a non-relativistic EFT. Note that at large box sizes, which 
are used in the study of the scattering processes, the characteristic center-of-mass momenta are small 
and hence working within the non-relativistic EFT is justified. 

An approach, which we use in the present paper, is closely related to the one of Refs. [12-16,20]. 
To be precise, we investigate the energy levels of the system with the quantum numbers of the A- 
resonance in a finite box with the size L. This is achieved by calculating the self-energy of the A in the 
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Figure 1: A schematic representation of the avoided level crossing in the presence of a narrow res- 
onance. The center-of-mass momentum of a two-particle pair is plotted against the size of a box L 
(arbitrary units). It is seen that in the vicinity of the resonance momentum pr a peculiar behavior of 
the energy levels is observed. The exact position of the resonance in this case could be easily pinned 
down by measuring the energy levels on the lattice. 

small scale expansion (SSE)# 5 at a finite volume and finding the poles of the propagator, which yields 
the parameterization of the spectrum of the Hamiltonian in a finite Euclidean box as a function of the 
variables M n and L, as well as the physical masses m/\(M n ), m^(M 7T ) and the irNA coupling constant, 
denoted g n NA- We further investigate the behavior of energy levels with respect to these parameters. 
The main question which is addressed here is whether the dependence of the energy levels on the 
above-mentioned parameters is pronounced enough in order to enable one to perform an accurate fit 
to determine the mass and width of the A-resonance. Here we would like to mention that, as it turns 
out later, the width of A is so large that the avoided level crossing is almost completely washed out 
from the energy spectrum. Therefore, the nice procedure shown in Fig. 1 does not immediately apply 
and, in order to be able to accurately extract the parameters of the A-resonance even in this case, 
one has to find an optimum fitting strategy. The discussion of this issue constitutes the main content 
of this work. 

The paper is organized as follows. In section 2 we consider the formalism that is used to calculate 
the self-energy of the A-resonance in the SSE, including a discussion of the infrared regularization 
procedure in a finite volume. Section 3 is devoted to the numerical calculations and the description 
of the fitting algorithm and the presentation of our results. In section 4 we introduce analytic param- 
eterizations of the energy spectrum at finite volume, which might be useful for performing fits to the 
lattice data. Finally, section 5 contains our conclusions. 

#5 The SSE is a phenomenological extension of chiral perturbation theory in which the nucleon-A mass splitting is 
counted as an additional small parameter. This quantity, however, does not vanish in the chiral limit. The framework 
of the SSE is laid out in detail in Ref. [30] 
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2 The formalism 



2.1 Extraction of the energy spectrum 

Let us consider the two-point function of two A-fields, given by Eq. (1), in a Euclidean box and 
choose the rest frame. After projecting out the spin-| state, this two-point function in momentum 
space develops poles at p™ = (iE n ,0). In coordinate space, the same two-point function at large 
Euclidean times t will be given by a sum of exponentials 

£<o|x M (x)x,(o)|o) L ->5>*V ex p(- £ w)- (2) 

x n 

where the subscript L indicates that the calculations are done at a finite volume L x L x L x Lf. We 
further assume that the size of the box in the time direction Lt is much larger than in the spatial 
direction L (eventually, L t — > oo) Consequently, studying the behavior of this two-point function at 
large Euclidean times, one is able to extract the energy levels E n = E n {L) of the positive parity, 
spin-| states, which are functions of the box size L. In case of a stable A-state, finite-size corrections 
to the lowest energy level vanish exponentially, E\ (L) — ttt-a = exp(— const • L) and in the large-L limit 
this level yields the value of the stable A-mass. This is, however, not the case for the decaying A, 
when the dependence of the energy levels on L is governed by a power rather than by an exponential 
law. The question is, whether one can extract the parameters of the A-resonance from the measured 
dependence of E n (L) on L. 

2.2 The Lagrangian and Feynman rules 

In order to parameterize the volume-dependent energy levels of the system in terms of the A-resonance 
parameters, we calculate the same two-point function in the SSE at finite volume. These calculations 
are similar to the recent study of the volume-dependent nucleon mass [31-35] - except for the fact that 
in the present paper we deal with an unstable particle. The calculations are performed by using the 
effective chiral Lagrangian, which explicitly contains pion, nucleon and A degrees of freedom [30, 36] 
and coincides with the Lagrangian used at infinite volume. The relevant terms (in Minkowski space) 
are displayed below: 

C = El (dHUdprf + X ] U + U^ X ) + • ■ ■ + C*N + ArA + £*na ; (3) 

c nN = 4^+4^+4^ + ... = ^A« +Ag+AS + ...}^v, 

A« = iP-m N+ ^ ^ 75 , 

AS = c 1 (x + }-—^({u^)D^D^ + h.c)+^(u^} + ---, 

Ag = <A (x + ) + <Ag + ...; (4) 

C^na = C^ NA + • • • , 4atA = g-KNA^aO^W^N + h.C. J (5) 
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ArA = 4i+^ + 41 + ■■■ = -^{(4i)^ + (41)^ + (^)^ + ■■■}0 , ' /, ^, 

( AS )^ = 0"° (^(<^« - ™a 4 /2 ) - ^ {7.7,3, (i Pa - m A 4 /2 )}) Of* + I ^ 755^ , 

4(mA) 2 ^ 



(A®)^ = ^%-{5 1 A Ao(x + ) + So A A^} + ---, (6) 

where (• • •) denotes the trace in flavor space. Throughout, we work in the isospin limit m u = md = m. 
The building blocks for the above Lagrangian are defined by 

U = u 2 , u fl = iu^df.Uu^ , = <9 M + ^ [u\d t 



9 <^«], 



X = 2B(s + ip) , x+ = u^X u ^ + , s = "il + 



1 



^ = - <t'u„> , D§ = Syr? - ie ijk (r k D») , u% = 8^ (7) 



and 



In these formulae, U = exp{ir • tt/F}, and ip^ represent the pion, nucleon and A interpolating 
fields, respectively. Further, rh^, % are the masses of the nucleon and the A in the chiral limit, 
and Ao =rn 1 \ — is the A-nucleon mass difference in this limit. The pion mass is given by 
M 2 = 2Bm + 0{rh 2 ). Furthermore, F, gA, flViVA denote the pion decay constant, the nucleon axial- 
vector constant and the 7riVA coupling constant in the chiral limit and q, Bf , a», B^ denote various 
low-energy constants (LECs). Finally, the projectors onto the isospin-| and isospin-^ subspaces are 
defined by 

The free A-propagator in d-dimensional space can be read off from the quadratic part of the Lagrangian 

l (0|T^ i (x)^(0)|0) = J ^ e-** S°,(p) f /2 ' lJ (10) 

with 



+ -A o {(^W + ^YVj ^4 — (^+^a)( j p 2 1 2 / V> (ii) 

v^T rh A < > (d-l)(m A ) 2 
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■ ^ a 1 ' b 0> b 1 

E L = ' = + < - + X 

A N 



Figure 2: The self-energy of the A at 0(e 3 ) in the SSE. 
where the projectors onto the spin-| and spin-^ states are defined as 

( P 21 2 )^ = ~^f=J~2 (P'P^v ~ P[xPv) , 
/pl/2v 

(Pi/V = ^ - V - V ■ ( 12 ) 
These projectors obey the following relations 

(P 3/2 )»x(P?/ 2 ) Xu = (PH 2 ),x(P 3/2 ) Xu = 0, 

(P 3/ \x(P 3/2 ) Xu = (P 3/ X> M'M = !' 2 - ( 13 ) 
2.3 The poles of the propagator 

The A propagator after inclusion of the self-energy diagrams (see Fig. 2) becomes 

d d p 



<0|2Vi(s)#(0)|0> = J — ^ e-*» S*(p) e 3/2 '^' (14) 



where 



SjlAp) = — (P 3/2 ),u + • • • , S L (p) =^Sf (p 2 ) + m A ^(p 2 ) , (15) 

"T-A - ^ + S L (p) 

and the ellipses stand for the terms with spin-^ projectors. The poles of the propagator are given by 
the zeros of the denominator 

m A - j*+ j*£fV)+m A £|V) =0- (16) 
We further introduce the quantities 

Sf(p 2 ) = £fV)-Re£,(p 2 ), * = 1,2, (17) 
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where =p"Ei(p 2 )+ ^(p 2 ) denotes the self-energy part of the A-resonance, calculated in the 
infinite volume (note that this quantity is complex below the decay threshold). Introducing now the 
physical A-mass (the real part of the pole mass) through the conventional definition (in the infinite 
volume) 

m A =m A (1 + ReSi(mi) + ReS 2 (m^)) + 0(p 4 ), (18) 

it can be seen that, at the order we are working, Eq. (16) can be rewritten by using the physical mass 
instead of the mass in the chiral limit 

m A - i>+ p't[{p 2 )+m A t%{p 2 )=0. (19) 

Performing the analytic continuation to Euclidian space p^ = (po>p) = (^P4>p) an d choosing the 
center-of-mass frame p^ = (a;, 0), from Eq. (19) we finally get 

m 2 A (l + t^-uj 2 )) 2 + lu 2 (1 - Ef 1 (-uj 2 )) 2 = . (20) 



2.4 Calculation of the self-energy at 0(e 3 ) 

Only pion-nucleon and pion-A loops, which are shown in Fig. 2, contribute to the quantities £f (— uj 2 ). 
The counterterm contribution is independent on L and thus drops out from the final expression. One 
obtains 



Sf {-uj 2 ) = £U-uj 2 ) + tf A (-uj 2 ) , i = 1, 2 , (21) 



where 



^(-uj 2 ) = [w 2 n {-uj 2 ) - % n {-uj 2 )} , t^ N (-uj 2 ) = ^ W? \-uj 2 ) (22) 

2£a("" 2 ) = ^{-fZ + (m 2 A+ u 2 )Wf(-u 2 ) + (m 2 A -^^ 

~ 3^ ( 3 < + u, 2 )# 2 A (-u, 2 ) - ^ (mi - c 2 )W 3 A (-u/ 



and 



t^ A (-uj 2 ) = ^[-f S + { m \+uj 2 )Wt{-uJ 2 )-2uj 2 Wt{-uJ 2 ) 



~ t^-^ + l^^-^}- (23) 

In the above expressions, the quantity Tq corresponds to the pion tadpole and W^{— uj 2 ), W^(—uj 2 ) 
to the ttN and ttA scalar loop functions, respectively (note that the tadpoles emerge from the one-loop 
diagrams shown in Fig. 2 after simplifying the numerators). An exact definition of these quantities 
is given below. Further, in these expressions one may take d = 4 because all ultraviolet divergences 
are contained in the infinite-volume integrals which, at the end, are included in the definition of the 
infinite- volume mass m A . Further, at the order we are working, the quantity F can be replaced by 
the pion decay constant F w . At the end, the self-energy part in a finite volume at 0(e 3 ) is expressed 
in terms of physical quantities only. 
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In the calculations at finite volume the loop integrals are replaced by infinite sums over discrete 
lattice momenta k n = 2im/L, n G Z 3 in spatial dimensions (see, e.g. [37]). To ease the notation, 
below we give the loop functions with the infinite-volume part included. In order to arrive at the 
quantities that enter Eqs. (22,23), one has first to isolate the infinite- volume part and then subtract 
it, e.g. 

TS = T% -TS^ . (24) 

Consequently, the expressions that are given below, are ultraviolet divergent and imply the application 
of some kind of regularization. Here we do not refer to a particular regularization procedure explicitly, 
since the divergent infinite-volume part will be always subtracted before the actual calculations are 
performed. 

The one-loop integrals that contribute to the self-energy of the A-resonance, are listed below 
(X = N,A). These are tadpole diagrams 

_fNM v l 



oo 27T L3 Z.. M 2 + k 2 



x _ fo- dk4 I 1 

J - / -KZTzl^ 



2vr L^m\+kr 

( x PaP A Ajn PaP/3 ™ f 00 dk 4 1 k a kp 

l^--^-r 2+ -^- Ta -- loo ^ ^MfTfcr 

( r PaPA M X , PaP/3 T X f°° dk ± 1 V- 

{ St #--p-) N * + ^ T * --J-^T^^TVJ (25) 



as well as the meson-baryon loop functions 
r°° dk 4 1 



thltXi 2\ _ f «fc4 1 \ ^ 1 

° l " ' ~ J-oo 2VT L 3 2. (M 2 + k 2 ){m 2 x + {f) _ kn?) > 

" W x ( — 2 "\ — [°° ^ 4 \ \ ^ ^£ 

Pa 1 ^ ] " loo ^ {Ml + kl)(m\ + (p- k n y) ' 

e PaPl3\ w X ( l ft,PaP/3 q X, , ,2x [°°dk 4 1 ^ 

6 ^~^n (_ " } + ~F~ 2 ( } " " 7-00 V M + + (P - ^n) 2 ) ' 



(M* + + <5^5 a ) Wf (-u, 2 ) + ^|f^ ^(-u; 2 ) 

/°° <i/c4 1 k a kpk a . . 

-oo ^TI^ ^ (M 2 + fc 2 )^ + _ fen) 2) > ^ ^ 

with the Euclidean 4-momentum (k n )^ = (&4,k n ). After subtracting the infinite-volume piece, the 
following recurrence relations between various functions can be obtained: 

m = \{Mlf--f-), N? = \{m\fx -f 2 x ), 

~ x , ^ + m x -Mt ~ x , TS f x 
W i^) ~ ^2 W o H" )-^2+^2' 

^ x (-^ 2 ) = Y^{ A( -^' m *' M2)# ° X( - w2) -^ 
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Wf(- 2 ) = ^ + % Ml W 2 X {-^) - JL { M ffi - n - miff + gj*} , (27) 
where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz + 2zx denotes the usual triangle function. 
2.5 Calculation of the scalar integrals 

The calculation of the tadpole graphs is straightforward and is carried out by using standard techniques 
(see, e.g. [33,34]). We demonstrate the method for the case of the pion tadpole. Using dimensional 
regularization to tame the ultraviolet divergence in Eq. (25), one gets 



° " L3 i-oc 27T ./ ^ k M 2 + k 2 + k 2 ^ d { 



27rn 

k 

L 



where the Poisson formula 



J_!l \ V e iLkj (281 

(2n) d M 2 + k 2 + k 2 Y ' 1 } 



+oo +oo 

S(x-n)= e2ninX ( 29 ) 

n=— oo n=— oo 

has been used to arrive at the second equality. Further, in this sum the term with j = corresponds 
to the infinite- volume integral. Separating this term, we finally obtain 



io - J (2n) d Ml + k\ + k 2 + ^ 4vr 2 Lj Jo ^ A 6 



To 



d d k 

M 2 ^ K^Lj) 

where j = |j| = + j 2 +i| and K u (z) denotes the modified Bessel function. 
The final result for the tadpoles is given by 

fn = M%y* KijMyLj) fX = mj L K x {m x Lj) 
47r 2 f^ M n Lj ' 4vr 2 ^ m x Lj ' 

~ M$ v K 2 (M*Lj) ~ x m x v K 2 (m x Lj) 

4n 2 ^ Q (M^Lj) 2 ' 2 47r 2 ^ (m x Lj) 2 ' { * > 

In the calculation of the meson-baryon loop functions one has to distinguish between two cases. In the 
7rA loop, the variable — lo 2 is below threshold. Using the Feynman parameterization, one may combine 
two denominators and then use the same technique as for the calculation of the tadpole contribution. 
As a result, one gets 



x, —uj 2 ' 



gA(x,-uJ 2 ) = (1 - x)M 2 + xm\ + x(l - x)uo 2 . (32) 
Note that for — J 2 close to m\ the function g&(x, -co 2 ) never vanishes in the integration region. 
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In contrast to the above example, the ttN loop can not be calculated by using the same method, 
because the variable —uj 2 can now be above the decay threshold A — > Ntt. To calculate this quantity, 
the following trick has been used. First, in order to avoid the ultraviolet divergence in the infinite sum 
over momenta, we have subtracted the integral at some scale uj 2 = fj 2 below threshold (one subtraction 
is enough for the convergence, but double subtraction enables one to achieve faster convergence). Now 
in the subtraction terms one is allowed to use the same technique as in the calculation of the irA loop, 
because — /U 2 is below threshold. This strategy is illustrated below in detail. The quantity Wq {— uj 2 ) 
which we are looking for is split into several terms 



W n (-uj 2 ) 
H^-uj 2 ) 
H 2 (-uj 2 ) 



W^(-u z )-ReW^(-u z ) 



L— >oo 



H^-uj 2 ) + H 2 (-uj 2 ) + H 3 (-uj 2 



[w n (-uj 2 ) - V) + (- 2 - 



Hsi-u 2 ) = -{Re<(-u; 2 )-<(V) + (- 2 -^)^<(-- 2 )L 2 _ 2 } 



L— >oo 



(33) 



The first term is the twice-subtracted infinite momentum sum, where the integration over is ex- 
plicitly performed 

1 ^ En + E w ; ! 

I? 



^i(- 2 ) = (- 2 -^E 



2E N E n uj 2 + (E N + E w ) 2 (p? + (E N + E n ) 2 ) 2 



E N = ^Jm 2 N + k 2 , E w = ^/M 2 + k 2 . 
The second expression corresponds to the subtraction term 

^2 I dxJ2(K (Ljy/ 9N {x, 



(34) 



H 2 {-u 2 



8vr 2 
(- 2 



2 x(l-x)Lj V ( T . 2T 

^V9N{x,-fi 2 ) \ v 



g N (x, -[i 2 



(1 - x)M 2 + xm 2 N + x(l - x)i? . 



(35) 



Furthermore, the remainder is included into the third term, which contains only quantities evaluated 
in the infinite volume 



H 3 (-uj 2 ) = - 



32it 2 uj 2 



In 



-uj 2 +m 2 N -M 2 + B t 



TT 



+ 



B„ 



16ttV 

uj 2 -11 2 
16ir 2 fi 2 

^arctan 



-uj 2 +m%- Ml - B^ 

,2 , _2 



+ In ■ 



-uj 2 — m 2 N + M 2 + B { 



—uj* — m 



N 



+ M 2 



1-bJ 



-fi 2 + m z N -M 2 -fi z -m N + M 2 

arctan ^ + arctan — 



B 



it 



B„ 



2\2 



1 + 



(uj 2 - n )(m N - Ml) ^ m 2 N f{m N + M 2 ) + (m 2 N - M 2 ) 



-li 2 + m 2 N 



B„ 



2uj 2 fi 2 
-Ml 

h arctan 



Ml 

-li 2 -m 2 N + Ml 



B 



B^ = A 1 / 2 (-w 2 ,m 2 v ,M 2 ), B„ = \V\- ll 2 ,m 2 N ,M%). 



(36) 
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2.6 Remarks 

a) In order to study baryon resonances, the volume should be taken much larger than in the case of 
stable particles since for excited states, the volume-dependent effects decrease only as powers of 
L and not by an exponential law as for stable particles. More precisely, the parameter L should 
be large enough, so that one could neglect all contributions of the type exp(— const • M^L) (see, 
e.g. [37,38]) as compared to the corrections that decrease according to the power law. 

b) The Lagrangians, given by Eqs. (3,4,5) and Eq. (6) do not contain the so-called "off-shell" 
pieces that do not contribute to observable quantities like the S'-matrix elements or transition 
currents. They could in principle contribute to the self-energy of the A, which is not an "on- 
shell" quantity. However, eliminating these terms at the diagrammatic level corresponds to 
canceling one of the propagators in the loop. The pertinent diagram turns into a tadpole, which 
vanishes exponentially. Therefore, for large values of L, where the exponential factors can be 
neglected, the off-shell couplings do not contribute to the self-energy. 

c) The same line of reasoning can be used to show that the contribution from the spin-i components 
of the A-propagator is irrelevant at finite volume. 

d) As seen from Eq. (32), the whole L-dependent part of the contribution from the 7rA loop is expo- 
nentially suppressed at large L. The same is true for all tadpoles. Even if we have retained them 
in the final expressions for completeness, the L-dependent part thereof can be safely neglected 
at any stage of the calculation. Indeed, we have checked numerically that the contribution from 
the 7rA loop is very small and does not affect the results. 

e) Covariant SSE calculations in the baryon sector at an infinite volume are performed by e.g. using 
infrared regularization (IR), which leads to a consistent power counting in the presence of the 
(large) baryon mass. This method has been also applied in the calculations at finite volume (see, 
e.g. [33]). The procedure, adopted in that work, amounts to merely extending the integration 
interval over the pertinent Feynman parameter from [0, 1] to [0, oo[, in a complete analogy to what 
is done at an infinite volume. In other words, in the context of the present problem it is equivalent 
to setting the N and A tadpoles to in Eq. (31) and to extending the integration over dx in 
Eq. (32) from to oo. It is, however, not immediately clear how this procedure can be generalized 
to the case when -co 2 is above threshold, see Eqs. (33,34,35,36). Our method is different from the 
one described above. Namely, we note first that applying IR in an infinite volume is equivalent to 
using ordinary dimensional regularization and changing the renormalization prescription. Next, 
going to a finite volume implies using the same Lagrangian, while replacing integrals by sums 
in the loops. In other words, our prescription reduces to using IR only in the infinite-volume 
self-energy whereas in the finite- volume piece S L (p) = S L (p) — S(p) ordinary dimensional 
regularization has been used. It is clear that in our case, unlike Ref. [33], the counterterms, 
corresponding to the above-mentioned additional renormalization, are L-independent. 

In the case of a stable particle changing the prescription from one to another amounts to intro- 
ducing corrections which vanish as exp(— const • mxL) with X = N, A at large L. This happens, 
in particular, for the volume-dependent nucleon mass, considered in Ref. [33]. It is clear that in 
this case the difference between the two prescriptions is physically irrelevant. It remains to be 
seen, what kind of statement can be made in case of an unstable particle. 
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Figure 3: The dependence of the energy levels in a finite box on the box size L for different values 
of the coupling constant g w NA- The avoided level crossing, which is clearly seen at small values of 
9-kNA (dashed lines), is washed out for the physical value of this coupling constant (solid lines). For 
comparison, we also display the free energy levels (dotted lines). It is seen that the energy levels in 
the presence of the interaction interpolate between different free energy levels. As expected, an abrupt 
change emerges in the vicinity of the resonance energy (the resonance position corresponds to the solid 
horizontal line). 



3 Numerical results 

Solving Eq. (20) for different values of L numerically we get the dependence of the energy levels E n (L) 
of the system, placed in a Euclidean box of size L. At this order, the energy levels depend on two 
parameters, namely the mass tba and the coupling constant g n NA- The term containing the coupling 
constant g± is exponentially suppressed and we assume that rn/v, M-^^F^ are determined independently 
in the same lattice calculations. We wish to investigate whether one can determine ttja and g n NA at 
a reasonable accuracy from a fit to the energy levels. 

The results of calculations of the energy levels for the physical values of all parameters are shown 
in Fig. 3 (solid lines). We use the following values for the particle masses: M n = 140 MeV, = 
940 MeV and rrij\ = 1210 MeV (the real part of the pole mass) and for the pion decay constant 
F n = 92.4 MeV. Further, we fix the physical value of the ttNA coupling constant as g^NA = 1-2 (at 
M w = 140 MeV). This corresponds to the value g\j^ = 13.5 GeV -2 {g-nNA = F^gA) obtained by 
using the pole approximation in dispersion relations, see Ref. [39]. Finally, we take g\ = 2.0 from 
Ref. [36] (nothing changes if we take g\ = 0). As we see, the coupling of the A to the 7r./V-system is 
so strong that the nice structure with the avoided level crossing has been almost completely washed 
out. It will resurface again, if the input value of g n NA is drastically reduced by hand (dashed lines 
in the same figure). This property, however, can not be meaningfully used in the fitting procedure. 
It is clear that, in order to perform the fit, another strategy, not linked to the identification of the 
resonance energy from the position of the avoided level crossing, should be looked for. 

In order to find such a strategy, we continue to study the structure of the energy levels, as well as 
the dependence on all available parameters. In particular, we start with varying ttja keeping all other 
parameters fixed. Consider for instance the lowest energy level shown in Fig. 4, where the dependence 
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Figure 4: The lowest energy level for different input values of toa and the physical values of g^NA 
and Mj. 

of the energy eigenvalue on L is plotted for different input value of m&. It is seen that the curves are 
almost linear in the variable L = M n L from the interval shown, with the tangent that monotonically 
decreases with decreasing tua- This property can be used for extracting uia- More precisely, one 
has to measure E\ on the lattice for different L and fit the curve, treating toa as a free parameter. 
The lattice data decide which of the curves in Fig. 6 has to be chosen. Since each curve corresponds 
to a particular value of the itia, the fit to the lowest level determines this parameter unambiguously 
(provided g n NA is known). 

Next, we consider the possibility of the determination of the irNA coupling constant. To this end, 
one has to find a quantity that could be maximally sensitive to g^NA and try to extract this coupling 
by a fitting procedure. Bearing in mind the level structure in the presence of a very narrow resonance 
(avoided level crossing), we expect that the difference of two energy levels Ei{V) — E\(V) can strongly 
depend on the width of the resonance. In Fig. 5 we plot the quantity Ei(V) — E\{V) against the 
variable L at a fixed value of ttia and varying the parameter g n NA- The avoided level crossing in this 
quantity - at small values of the parameter g-nNA ~ is seen as a sharp minimum near the value of L 
where the crossing takes place and the value of the function at the minimum determines the width. 
As evident from Fig. 5, even at the physical value of g n NA one may observe a remnant of the avoided 
level crossing - a plateau, which disappears if g-wNA increases further. Given a rather pronounced 
dependence of the quantities plotted in Fig. 5 on the input value of g n NA, one may expect that fitting 
would allow one to determine this coupling constant with a reasonable accuracy. 

We are now in a position to describe our proposal for determining the parameters of the A- 
resonance tba and g n NA from the lattice data, which is expected to work, even if the width of the 
resonance is not very small. In brief, we propose to fit the first few energy levels, measured on the 
lattice, to the calculated energy levels, which are parameterized by the free parameters itia and g n NA- 
Further, it could be advantageous to carry out this procedure iteratively. Fix first the decay constant 
to some input value and determine the A-mass from the L-dependence of the lowest energy level, see 
Fig. 4. With the newly determined A-mass plot the difference E2(L) — E\{L) and fit the parameter 
g^NA to the data, see Fig. 5. Repeat the procedure until convergence is achieved. 
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Figure 5: The difference of the first two energy levels for different values of g n NA and the physical 
value of M n . For small values of the coupling constant there is a dip, corresponding to the avoided 
level crossing. A plateau is clearly visible even for the physical value of g n NA- 

Finally, we wish to comment on the dependence of the energy levels on the other parameter which 
is at our disposal, namely the quark (pion) mass. Since this structure depends only on the physical 
nucleon and A masses, in the present (exploratory) study of the problem we have restricted ourselves 
to the 0(e 2 ) expressions of the baryon masses 

m N = mfy* - 4ci (M2 - (Mf^) 2 ) , m A = mf - 4 ai (M 2 - (M^ s ) 2 ) (37) 

and take c\ = —0.9 GeV -1 , a\ = —0.3 GeV -1 [36] (note that c\ and a\ enter only through the nucleon 
and A masses in the infinite volume, so the choice of particular numerical values for these constants 
does not affect our conclusions). Moreover, we neglect the pion mass dependence of the coupling 
constant F^, since the pertinent correction arises at higher order. The above simple case perfectly 
models the real situation: the mass difference between to a and + M n monotonically decreases 
with the increase of M v and vanishes at around M n = 210 MeV. Below threshold, the structure 
of the energy levels and the dependence on the parameters mA, g n NA is similar to the case with 
M v = 140 MeV. After crossing the threshold from below, the A is stable and one expects that the 
finite-volume corrections to the lowest energy level get exponentially suppressed. 

Fig. 6 clearly illustrates this pattern. In this figure, the lowest-order energy level is plotted at 
three different values of the pion mass and the physical value of the decay constant. As we see, the 
curve is rather smooth in all cases and monotonically flattens as the decay threshold is approached 
from below. This property can be discussed in a more quantitative fashion. Namely, fitting the level 
energy in Fig. 6 E\ = Ei/rriN to the variable L within the interval 3 < L < 8.3 with a linear function 
Ei = A + BL, for the different values of the mass gap u>o = m& — — M n we get: (A,B) = 
(1.34,-0.019) for lj = 130 MeV, (A, B) = (1.30,-0.009) for uj = 78 MeV, (A, B) = (1.25,-0.006) 
for ujq = 20 MeV. It is seen that, while the parameter A remains almost stable, the linear coefficient 
decreases monotonically. It is important to observe that there is no sign of irregularity in the quantity 
Ei(L), when one crosses the threshold. Finally, we would like to mention that here the distance at 
which the chiral extrapolation takes place is not bound from below and could be taken smaller than 
the distance to the decay threshold. Consequently, the extrapolation error could be reduced. 
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Figure 6: The lowest energy level for different values of M n and the physical value of g n NA- The 
curves correspond to the following values of ujq = m& — tun — M^: a) solid line: uq = 130 MeV, b) 
dashed line: uj = 78 MeV, c) dotted line: uj = 20 MeV. 



4 Analytic parameterization 

The dependence of the energy levels E n {L) on the parameters ro^, g n NA an d L is given by the 
numerical solution of Eq. (20). If one implements the above fitting procedure in practice, it would 
be useful to have a simplified algebraic expression for the energy levels, where the dependence on the 
parameters is explicit. Below we demonstrate how such an expression can be derived. Starting from 
Eqs. (21,22,23), we neglect all contributions that are exponentially suppressed in L. Namely, the irA 
loop is neglected altogether as well as various tadpoles. The functions £j(— uj 2 ), which appear in the 
self-energy part, are then proportional to the loop function Wq {— uj 2 ), defined in Eq. (26). Then, 
Eq. (20) simplifies to 
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(38) 



Further, since the A is a P-wave state, the lowest singularity in the self-energy corresponds to the 
contribution of the term with n 2 = 1 (it can be easily checked that the singularity at n 2 = cancels 



with the factor A(- 
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, M 2 ) in the numerator). Isolating the singularity at n = 1 in the function 
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and the function Rq (— uj ) is regular in the vicinity of uj 
rewritten as 
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Further, Eq. (20) can be 



(41) 



16 




ML ML 

7T 7T 



Figure 7: Approximate solution for the first two energy levels, Eq. (43), for physical (left panel) and 
reduced (right panel) values of the coupling constant g n NA- The approximate solution qualitatively 
reproduces the features of the exact solution. Namely, the avoided level crossing is clearly visible in 
the right panel, as well as an almost stable A-state at M n = 200 MeV (left and right panels). 

where 

S < L) = 16^5? { (E<11 +m " )2 - ^ 2 } A « £(1) > 2 ' ro »> M -> Z^Ii« (42) 

and the function r(y/ —uj 2 ) is reg ular at V~u 2 = 

If now we use the approximation r(y— uj 2 ) = 0, we arrive at a quadratic equation, whose solution 

Ei, 2 (L) = ~ |m A + ± y/(m A - £«) 2 + 4<?(L)j (43) 

gives the position of the first two energy levels. In Fig. 7 these levels are plotted at the physical value 
of Q-kNA j as well as for the case of a smaller value. It is amusing that such a simple solution reproduces 
the gross features of the exact solution quite well. For example, the avoided level crossing is clearly 
visible at smaller values of g^NA- Finally, we note that in Ref. [16] a very similar equation is derived 
in a simple two-channel quantum-mechanical model. 

The approximate solution, given by Eq. (43), contains an explicit dependence on the parameters 
L, raA and g v NA- For this reason, it is easier to use it for performing the fit to the lowest energy 
levels. Of course, from the quantitative point of view, the quality of the approximation is still not 
satisfactory. It is however obvious that the accuracy can be systematically improved by including the 
next nearby singularities as well as regular contributions. We do not display the pertinent formulae 
here. It is clear that, at the end, the accuracy of the analytic parameterization, which is used to 
analyze the real lattice data, should be matched to the precision of this data. 

5 Conclusions 

a) In this paper we present the results of calculations of the pole structure of the correlator of 
two A-fields in a finite Euclidean box. It has been argued that, calculating the first few energy 
levels E n (L) in terms of the resonance parameters toa and g v NA within the SSE and fitting the 
lattice data at finite volume, an extraction of these parameters at a reasonable accuracy may be 
possible. This statement constitutes the main result of the present paper. 
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b) The main question that remains is, how the result will be affected by higher-order corrections in 
the chiral expansion. It can be for instance shown that, parameterizing pion-nucleon scattering 
matrix by the pure s-channel A-pole term, Eq. (20) can be rewritten in a form similar to 
the Liischer's master formula, which expresses the displacement of an energy level through the 
scattering phase shift. The higher-order corrections in the SSE contribute to the non-resonant 
background and are expected to be moderate. Of course, such heuristic arguments can not be 
a real substitute of explicit calculations. Note that such calculations at 0(e 4 ) are already in 
progress and the results will be reported elsewhere, including a detailed error analysis in SSE at 
this order [40] . One expects that these calculations shed light on the question of the convergence 
of the chiral expansion. 

c) Our approach is closely related to the method proposed originally by Liischer and developed in 
number of subsequent publications [12-16,25]. Note, however, that the Liischer formula, relating 
the energy levels of a system in a finite volume to the scattering phase shifts, is valid in general 
beyond the chiral expansion, thus avoiding the above-mentioned problem of convergence. On 
the other hand, our result contains an explicit parameterization of the energy levels in terms 
of roA, QnNA and can be used, in addition, to study the quark mass dependence of the energy 
levels. This is important because the first lattice data which will appear below ttN threshold, 
will probably still correspond to the pion mass higher than the physical value. 

We now note that, due to the condition M W L » 1, the characteristic 3-momenta of the system 
p <C Mjr and the non-relativistic approach must be applicable - the processes with a mass gap 
~ M n or higher are suppressed exponentially. One therefore expects that, for a sufficiently large 
L these two approaches overlap and the results are complementary to each other. At present, 
we are investigating the problem in detail within non-relativistic EFT, aiming to explicitly 
demonstrate this relationship that, in turn, will enable one to choose an optimum strategy for 
determining the resonance parameters from the lattice data in the future [41]. 
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